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ABSTRACT 

Joule  heating  in  the  upper  atmosphere  is  the  most  important  energy  dissipation  process  between  the 
magnetosphere  and  the  ionosphere.  In  this  paper  we  examine  the  various  terms  in  the  equation  that  governs 
Joule  heating. 

(1)  The  ionospheric  electric  field  seasonal  dependence  is  examined.  It  is  found  that  the  shape  of  the 
ionospheric  plasma  convection  cells,  and  the  latitude  of  the  reversal  from  sunward  to  antisunward 
convection,  are  seasonally  dependent.  Statistical  averages  of  the  square  of  the  ion  velocity  show  a  maximum 
in  fell,  and  a  minimum  in  summer. 

(2)  Pedersen  conductivities  at  F  region  altitudes  are  examined  using  Chatanika  and  Sondrestrom  radar 
data.  It  is  shown  that  during  solar  minimum  conditions,  the  F  region  contributes  less  than  20%  to  the  total 
height  integrated  Pedersen  conductivity  Ep.  In  contrast,  during  solar  maximum  conditions  the  contribution 
to  E,,  from  solar  produced  F-region  ionization  can  be  60%. 

(3)  The  importance  of  the  neutral  wind  term  in  Joule  heating  calculations  is  illustrated  using  a  specific 
example.  The  Joule  heating  calculated  by  including  the  neutral  wind  term  is  2  to  4  times  smaller  than  that 
calculated  without  the  neutral  wind.  But  the  reverse  can  also  be  true,  as  shown  during  a  period  when  the 
neutral  wind  played  the  role  of  a  dynamo  in  the  ionosphere/  magnetosphere  current. 


INTRODUCTION 

It  is  generally  accepted  that  the  solar  wind  energy  that  reaches  the  magnetosphere  is 
dissipated  by  three  main  processes:  (1)  injection  of  particles  into  the  ring  current  and  further 
dissipation  by  charge  exchange,  (2)  Joule  heating  in  the  ionosphere  and  thermosphere,  and  (3) 
particle  precipitation.  The  energy  transfer  rate  into  the  ring  current  is  often  larger  than  that  of  the 
other  two  processes.  At  Units,  however,  the  Joule  heating  rate  can  be  the  largest  [Baumjohann  and 
Kamide,  1984],  Although  particle  precipitation  in  the  auroral  oval  may  dissipate  energy  at  a  rate 
comparable  to  Joule  dissipation  [Evans  et  al.,  1977;  Vickrey  et  al.,  1982],  globally,  the  dissipation 
by  Joule  heating  is  significantly  more  important  than  by  particle  precipitation  [Baumjohann  and 
Kamide,  1984].  The  effect  on  the  neutral  atmosphere  of  the  Joule  dissipation  has  far-reaching 
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consequences.  Joule  dissipation  heats  the  neutral  gas  and  causes  upwelling  of  atomic  oxygen 
which  is  then  transported  to  lower  latitudes.  This  process  changes  the  global  configuration  of 
thermospheric  neutral  winds  [Roble  et  al.,  1977;  Prftlss,  1980].  As  a  result,  the  Joule  heating  rate 
is  one  of  the  most  fundamental  parameters  needed  in  studies  of  the  coupling  mechanisms  between 
the  solar  wind  and  the  thermosphere-ionosphere- magnetosphere  system. 


The  height-integrated  Joule  heating  rate,  in  the  reference  frame  of  the  neutrals,  is  governed 
by  the  equation 


(h)(E  +  U  X  B)2dh 


(1) 


where 


op  is  the  Pedersen  conductivity 

E  is  the  electric  field  calculated  from  the  ionospheric  plasma  drift  V  by  E  =  -  V  X  B 
U  is  the  neutral  wind 


B  is  the  magnetic  field 

h\  and  hi,  the  limits  of  integration,  encompass  the  E  and  F  region  altitudes. 

Incoherent  scatter  radars  are  the  only  ground-based  instruments  capable  of  measuring  all  the 
quantities  in  Equation  [1],  They  can  easily  measure  the  electron  density,  electric  field  and  neutral 
wind  vectors.  They  can  also  measure,  although  with  more  difficulty,  the  ion  composition  and  ion- 
neutral  collision  frequency  (Lathuill&re  et  al.,  1983;  Johnson  and  Wickwar,  1987). 

Because  an  excellent  review  paper  has  recently  been  published  on  energetics  of  the 
thermosphere  (Killeen,  1987),  we  do  not  review  the  subject  of  Joule  heating  again  in  this  paper. 
Instead,  we  concentrate  on  some  specific  aspects  of  each  term  in  Equation  [1].  We  first  consider  the 
electric  field,  and  show  that  seasonal  change  is  significant  in  the  electric  field  pattern  and 
intensity  observed  from  Sondrestrom  (de  la  Beatyardtere  et  al.,  1988).  We  then  consider  the 
conductivities,  and  show  that  during  solar  maximum  the  solar-produced  F  region  ionization  can 
contribute  more  than  50%  of  the  total  height-integrated  Pedersen  conductivity.  Finally,  we  make 
some  remarks  concerning  the  importance  of  the  neutral  wind  contribution  to  Joule  heating. 

In  this  paper,  data  from  the  Sondrestrom  and  Chatanika  incoherent  scatter  radars  (Kelly, 
1983,  Leadabrand  et  al.,  1972)  are  used.  At  75°  invariant  latitude,  Sondrestrom  is  usually  in  the 
polar  cap  during  the  night.  It  is  in  the  cleft,  cusp  and  auroral  zone  during  the  day.  At  65°  invariant 
latitude,  Chatanika  is  usually  in  the  auroral  zone  during  the  night,  and  equatorward  of  the 
auroral  oval  during  the  day. 


SEASONAL  DEPENDENCE  OF  HIGH-LATITUDE  IONOSPHERIC  CONVECTION 


The  dependence  of  high-latitude  electric  field  on  geomagnetic  activity  and  IMF  configuration 
has  been  studied  for  many  years,  using  data  from  satellites  and  ground-based  instruments 
(Heppner,  1977;  Holt  et  al.,  1987;  Heppner  and  Maynard,  1987;  Foster  et  al.,  1971,  1986;  Friis- 
Christensen  et  al.,  1985;  de  la  Beapjardifcre  et  al.,  1985,  1986;  de  la  Beaiyardi&re  and  Wickwar, 
1986;  Holt  et  al.,  1987).  Little  is  known,  however,  about  possible  seasonal  variations  of  electric 
fields.  This  is  surprising,  because  it  has  been  known  for  decades  that  the  various  indices  of 
magnetic  activity  are  seasonally  dependent,  with  maxima  in  spring  and  fall  (Corlie,  1912;  Russell 
and  McPherron,  1973,  and  references  therein).  Similarly,  high-latitude  ionospheric  currents,  and 
field  aligned  currents  have  maxima  in  summer  (Matsushita  and  Xu,  1982;  Campbell,  1982;  Fujii 
et  al.,  1981).  To  our  knowledge,  a  systematic  analysis  of  the  seasonal  dependence  of  ionospheric 
electric  fields  has  never  been  done  before  the  work  described  here  and  in  de  la  Beaqjardidre  et  al. 


(1988).  Seasonal  variations  in  plasma  convection  are  important  to  assess  in  the  context  of  Joule  lty  eodoa 
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heating  calculations  since  the  Joule  heating  rate  is  approximately  proportional  to  the  square  of 
the  ion  velocity. 

We  have  binned  and  averaged  all  the  Sondrestrom  data  taken  during  experiments  lasting 
more  than  7  hours.  The  observations  spanned  five  years,  roughly  centered  around  solar  minimum, 
from  March  1983  to  July  1988.  More  than  130  days  were  averaged.  The  most  salient  points 
revealed  by  these  studies  and  illustrated  by  Figures  1  to  3  can  be  summarized  as  follows: 

•  The  equinox  and  summer  convection  velocities  show  a  smooth  two-cell  pattern,  whereas  the 
winter  convection  velocity  is  comparatively  disorganized. 

•  The  winter  East- West  (E-W)  component  of  the  velocity  is  the  largest  in  the  dawn  and  noon 
sectors. 

•  The  summer  North-South  (N-S)  component  is  the  smallest  in  the  noon  and  dusk  sectors. 

•  There  is  a  phase  difference  in  the  daily  variations  of  the  E-W  velocity  component,  the  winter 
months  leading  equinoctial  and  summer  months.  Especially  on  the  day  side,  this  phase 
difference  can  be  interpreted  as  a  clockwise  rotation  of  the  convection  pattern. 

•  In  the  noon  portion  of  the  afternoon  cell,  the  latitude  of  the  convection  reversal  from  sunward 
to  antisunward  is  highest  during  the  summer  months. 

•  In  the  dawn  convection  cell  the  latitude  of  the  convection  reversal  from  sunward  to 
antisunward  is  highest  during  the  winter  months. 

We  now  discuss  some  of  these  points  in  more  details.  Figures  la  through  Id  show  the  average 
velocities.  Thick  and  thin  vectors  are  used  to  distinguish  west  and  east  velocities,  respectively.  The 
bin  size  is  1°  by  1  hour,  covering  the  invariant  latitude  interval  from  67°  to  82°.  The  number  of 
points  in  each  bin  is  typically  200  to  300.  Comparing  summer  and  winter  data  (Figures  lb  and  d), 
we  see  that  the  average  convection  is  somewhat  smoother  and  more  regular  in  summer.  This  is  not 
unexpected  because,  on  individual  winter  days,  and  especially  on  the  nightside,  the  convection  has 
previously  been  reported  to  be  highly  irregular  (Heppner,  1977).  This  is  evidenced  in  the  winter 
convection  pattern  shown,  for  example,  by  Heelis  and  Hanson  (1980)  and  de  la  Beaujardidre  et  al. 
(1985).  However,  the  averaging  of  winter  data  presented  here  smoothes  out  much  of  the 
irregularity,  and  the  winter  drifts  appear  to  maintain  a  somewhat  erratic  but  fairly  well  organized 
two-cell  convection  pattern. 

Figures  2a  and  2b  show  the  individual  components  of  the  plasma  drift  as  a  function  of  time  for 
all  seasons.  (The  data  for  this  figure  and  the  following  one  correspond  to  a  2°  by  2  hour  grid  size.) 
The  winter  E-W  component  is  largest  between  05  and  15  LT,  i.e.,  in  the  dawn  and  noon  local  time 
sectors.  The  summer  N-S  component  is  smallest  between  09  and  17  LT,  i.e.,  in  the  noon  and  most  of 
the  dusk  local  time  sectors. 

The  latitude  at  which  the  convection  reverses  from  sunward  to  antisunward  is  displayed  in 
Figure  3.  For  this  figure,  the  2°  X  2  hr  averaged  data  were  used.  The  latitude  of  this  reversal  has 
an  opposite  dependence  on  season  in  the  dusk  and  dawn  cells.  In  winter,  the  dusk  cell  has  its 
convection  reversal  at,  or  equatorward  of,  the  other  seasons.  In  summer,  the  dawn  cell  has  its 
convection  reversal  at,  or  equatorward  of,  the  other  seasons.  In  other  words,  the  dusk  convection 
reversal  is  at  lowest  latitudes  during  the  winter,  while  the  dawn  convection  reversal  is  at  lowest 
latitudes  during  the  summer.  These  latitudinal  differences  are  most  pronounced  in  the  noon  and 
midnight  sectors. 

Note  that,  for  the  summer  curves,  the  afternoon  reversal  is  beyond  the  field  of  view  at  noon 
when  it  is  poleward  of  81°.  At  midnight,  the  reversal  is  also  beyond  the  field  of  view,  but  it  is 
equatorward  of  68°.  The  winter  reversal  is  located  equatorward  of  the  summer  reversal  in  the 
afternoon  cell  and  poleward  in  the  morning  cell. 

Because  Joule  heating  is  approximately  proportional  to  the  square  of  the  electric  field,  or 
equivalently,  to  V*,  the  individual  drifts  were  also  squared  and  averaged  according  to  season.  The 
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results  of  averaging  all  latitudes  together  are  displayed  in  Figure  4  as  a  function  of  local  time.  The 
daily  variations  show  three  maxima  which  occur  (1)  one  to  two  hours  before  local  midnight;  (2)  at 
dawn;  and  (3)  one  to  two  hours  after  local  noon.  (At  Sondrestrom,  magnetic  local  time  precedes 
local  time  by  about  one  hour.)  Except  during  summer,  the  largest  values  are  those  that  occur 
before  noon  and  before  midnight.  Part  of  these  variations,  however,  are  due  to  the  offset  between 
the  geomagnetic  and  geographic  poles.  Over  the  course  of  one  day,  the  Sondrestrom  radar  samples 
differing  parts  of  the  convection. 

As  stated  by  Foster  et  al.  (1983),  the  seasonal  differences  in  the  square  of  the  velocity  are 
small.  However,  close  examination  of  Figure  4  reveals  that  the  small  difference  between  seasons  is 
probably  significant.  By  averaging  over  all  local  times  we  find  that  V2  is  the  smallest  in  summer 
(V2  =  1.6  X  106  m2  S'2)  and  largest  in  fall  (V2  =  1.8  X  105  m2  s-2).  Spring  and  winter  are  about 
equal.  In  fact,  the  same  seasonal  dependence  is  also  apparent  in  Foster  et  al.  (1983),  who  binned 
the  square  of  the  electric  field  measured  from  the  Atmosphere  Explorer  C  (AE-C)  satellite  during 
two  time  intervals  close  to  the  previous  solar  minimum.  By  inspecting  their  figure,  reproduced 
here  as  Figure  5,  it  is  clear  that  the  squared  electric  field  (E2)  is  largest  in  fall  and  smallest  in 
summer.  Also,  this  difference  is  more  pronounced  during  the  quiet  periods  (Kp  0  to  3)  than  during 
the  active  periods  (Kp  3  to  6).  And  indeed,  these  authors  state  that,  overall,  the  E2  is  20%  lower  in 
summer  than  in  winter  and  equinox. 

The  differences  in  E2  have  repercussions  in  terms  of  Joule  heating.  Neglecting  the  neutral 
wind  effects,  Foster  et  al.  (1983,  Figure  7)  find  that  during  quiet  times,  the  Joule  heating  rate 
integrated  over  all  latitudes  is  maximum  in  fall.  The  higher  conductivities  resulting  from  the 
larger  solar  illumination  during  summer  are  not  sufficient  to  offset  the  fact  that  E2  is  smaller  in 
summer  than  in  the  fall. 


CONDUCTIVITIES  FROM  SOLAR-PRODUCED  IONIZATION 
The  Hall  and  Pedersen  conductivities  are  given  by: 
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(3) 


Ne  is  the  electron  density 

e  is  the  charge  of  electron 

yen  is  the  electron-neutral  collision  frequency 

v0,  \0y  Vf,0  are  the  ion-neutral  collision  frequencies 

Qt  is  the  electron  gyrofrequency 

a0,  a0j,  am  are  the  ion  gyrofrequencies 

Po>Poj>Pn2* re  the  relative  ion  concentrations 

Recently,  Brekke  and  Hall  (1988)  reviewed  conductivity  measurements  and  empirical  models. 
Listing  the  various  terms  in  Equations  [2]  and  [31,  they  described  the  simplifying  assumptions 
that  are  often  made  in  calculating  the  Hall  and  Pedersen  conductivities.  In  particular,  these 
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authors  stressed  that  the  collision-frequency  models  adopted  are  of  great  importance  when 
deriving  the  conductivity. 

Brekke  and  Hall  (1988)  also  helped  clarify  the  dependence  of  conductivities  on  the  solar  zenith 
angle  (x).  While  several  authors  have  stated  that  the  Pedersen  and  Hall  conductivities  (Ep  and  Eh) 
are  proportional  to  cosx,  others  have  argued  that  conductivities  are  proportional  to  (cosx)0-5,  or, 
more  generally,  to  (cosx)n.  Brekke  and  Hall  (1988)  showed  that  the  conductivities  are  better 
described  using  an  expression  where  the  two  terms  in  cosx  are  added.  They  obtained  for  the 
Pedersen  and  Hall  conductivities  an  expression  of  the  form:  E  =  a  cosx  +  b(cosx)0-5-  Adding  these 
two  terms  is  justified  by  the  fact  that,  in  the  lower  E  region,  the  ionization  recombination  process 
is  dominated  by  dissociative  recombination,  so  the  electron  density  is  proportional  to  cosx;  whereas 
in  the  upper  E  region  and  in  the  F  region,  charge  exchange  processes  dominate  and  the  electron 
density  is  proportional  to  (cosx)0-5. 

In  another  recent  article,  Rasmussen  et  al.  (1988)  studied,  among  other  things,  the  height- 
dependent  Pedersen  conductivities.  They  calculated  how  much  of  the  height-integrated  Pedersen 
conductivity  originates  from  F  region  altitudes  (above  —170  km)  and  how  much  originates  from  the 
E  region.  These  authors  showed  that,  during  solar  maximum  conditions,  the  F  region  can 
contribute  more  than  previously  assumed.  They  derived  the  values  of  the  ionospheric  densities  and 
conductivities  from  the  combination  of  a  photochemical  equilibrium  model  at  E-region  altitudes, 
and  the  full  F-region  model  including  transport  and  diffusion  at  higher  altitudes.  Under  solar 
maximum  conditions,  they  concluded,  the  contribution  to  Ep  from  the  ionosphere  above  170  km 
can  reach  40%  during  day-time  and  80%  at  dusk.  Comparison  of  their  ionospheric  density  model 
with  incoherent  scatter  data  shows  that  in  the  altitude  range  110  to  200  km,  the  predicted  values 
are  underestimated  by  40  to  50%.  It  is  therefore  important  to  assess,  using  actual  data,  the  F 
region  contribution  to  height  integrated  conductivities. 

We  have  examined  Chatanika  and  Sondrestrom  ionospheric  density  measurements  and 
calculated  the  conductivities  (Equations  [2]  and  [3])  without  making  the  approximations  that  have 
often  been  made  in  the  past.  The  procedure  was  the  same  as  used  for  the  Chatanika  data  in 
Rasmussen  et  al.  (1988),  and  we  looked  at  all  the  days  that  were  analyzed  by  these  authors.  All  the 
collision  terms  between  electrons,  NO+ ,  C>2+ ,  0+  and  N2,  O2,  O  were  included,  using  the  collision 
frequencies  in  Schunk  and  Nagy  (1980).  The  neutral  atmosphere  was  obtained  from  the  MSIS-83 
model  (Hedin  et  al.,  1983).  The  results  are  shown  in  Figures  6a  through  6c  where  the  solid  and 
dashed  lines  represent  Pedersen  conductivity  integrated  up  to  500  and  180  km  altitudes, 
respectively.  Figure  6a  is  for  solar  cycle  minimum  (Sio.7=  69);  Figures  6b  and  6c  are  for  solar  cycle 
maximum  (S10.7  =  197  and  180,  respectively). 

Qualitatively,  these  calculations  agree  with  Rasmussen  et  al.  (1988).  During  solar  cycle 
minimum  the  difference  between  the  two  curves  is  very  small,  indicating  that  the  F-region 
contribution  to  Ep  is  small,  less  than  —10%.  During  solar  cycle  maximum,  the  difference  can  be 
larger.  Figure  6b  shows  a  large  F-region  contribution,  reaching  60%,  during  the  daytime  hours. 
Although  it  is  daytime,  the  sun  has  a  razing  incidence  and  the  solar  zenith  angle  always  exceeds 
80°  with  the  result  that  the  rate  of  EUV  photoionization  is  small  in  the  E  region  and  large  in  the  F 
region.  This  corresponds  to  the  dusk  conditions  described  by  Rasmussen  et  al.  (1988)  when  the  F- 
region  contribution  is  largest.  However,  this  contribution  is  smaller  than  that  obtained  from  the 
model  calculation,  which  was  80%.  Figure  6c  shows  a  much  smaller  F-region  contribution  than 
Figure  6b,  about  20%.  In  this  case  the  solar  zenith  angle  is  smaller,  giving  rise  to  a  larger  E-region 
contribution  to  the  Pedersen  conductivity.  This  is  again  in  qualitative  agreement  with  Rasmussen 
et  al.  (1988),  but  as  in  Figure  6b,  the  F-region  contribution  is  smaller  than  that  obtained  from  the 
model  calculation. 

In  order  to  asses  the  relative  importance  of  each  element  in  the  conductivity  calculation 
(Equations  [2]  and  [3]),  we  show  in  Figures  7a,  7b,  and  8  the  height  profiles  of  the  Hall  and 
Pedersen  conductivities  and  of  the  electron  density.  The  time  is  local  noon,  and  the  day  selected  is 
February  4,  1981,  the  day  illustrated  in  Figure  6b  for  which  the  F  region  contribution  to  the 
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Pedersen  conductivity  was  large.  Figure  7a  shows  that  the  contribution  to  the  Pedersen 
conductivity  from  the  electrons  (the  first  term  in  Equation  [2])  is  negligible  at  all  altitudes  above 
90  km.  The  contribution  from  the  molecular  ions  (last  two  terms  in  Equation  [2])  dominates  up  to 
about  190  km.  The  contribution  from  atomic  oxygen  ions  (second  term  in  Equation  [2])  dominates 
above  190  km.  This  is  because  the  electron  density  was  extremely  large  in  the  F  region  (1.8  X  106 
el/cm3),  whereas  it  was  very  small  in  the  E  region  (5.3  X  104  el/cm3). 

The  components  of  the  Hall  conductivity  are  displayed  in  Figure  7b.  The  Hall  conductivity  is 
the  difference  between  the  contribution  from  the  electrons  (first  term  in  Equation  [3])  and  from  the 
ions  (three  other  terms  in  Equation  [31).  The  conductivity  peaks  at  112  km,  slightly  below  the  E- 
region  peak  at  120  km.  Below  the  peak,  oh  depends  mostly  on  the  electron  term;  above  the  peak, 
molecular  ions  determine  the  shape  of  the  curve.  The  contribution  from  atomic  oxygen  is 
negligible  at  all  heights. 


NEUTRAL  WIND 

The  final  parameter  we  would  like  to  comment  on  is  the  neutral  wind,  U.  Equation  1  shows 
that  the  Joule  heating  of  the  neutral  atmosphere  is  proportional  to  the  square  of  the  vector 
difference  between  the  ion  drift  and  the  neutral  wind.  Because  few  techniques  can  provide  the  V 
and  U  vectors  simultaneously,  and  because  U  is  generally  smaller  than  the  E  X  B  drift,  the 
neutral  wind  is  often  ignored.  However,  its  contribution  can  be  quite  important.  In  the  F  region, 
ion  drag  tends  to  diminish  the  size  of  V-  U;  and  in  the  E  region,  the  wind  can  be  very  large  during 
geomagnetically  active  conditions  (Brekke  and  Rino,  1978;  Johnson  and  Wickwar,  1987;  Baron 
and  Wand,  1983).  Depending  on  the  magnitudes  of  the  ion  and  the  neutral  velocities,  as  well  as  the 
time  history  of  a  particular  event,  the  Joule  heating  rate  can  be  either  reduced  or  enhanced 
relative  to  the  U  =  0  value. 

The  incoherent  scatter  technique  provides  the  electric  field  at  all  heights.  This  technique  also 
allows  both  components  of  U  to  be  determined  in  the  E  region,  and  the  meridional  component  of  U 
to  be  determined  in  the  F  region  (Johnson,  1988).  Joule  heating  rates  derived  from  Chatanika 
radar  measurements  using  Equation  [1]  are  displayed  in  Figure  9.  These  observations  show  an 
example  of  the  magnitude  of  the  neutral  wind  effect  on  a  highly  disturbed  day  characterized  by 
large  electric  fields  and  intense  particle  precipitation.  The  two  curves  on  this  figure  show  the  local 
Joule  heating  rates  at  115  km  altitude.  The  solid  line  is  derived  from  Equation  1,  and  the  dashed 
line  computed  by  setting  the  neutral  wind  equal  to  zero.  The  observations  correspond  to  the  period 
illustrated  in  Figure  6c  (June  10-11, 1980).  As  this  particular  example  shows,  the  Joule  heating 
rate  that  includes  the  neutral  wind  effect  is  generally  2  to  4  times  less  than  the  rate  derived  by 
assuming  U  =  0.  The  decrease  is  highly  variable,  however,  and  at  times  approaches  or  exceeds  a 
factor  of  10.  After  20  UT  (10  LT)  the  Joule  heating  rate  derived  by  assuming  U  =  0  is  actually 
exceeded  by  that  obtained  by  including  the  neutral  wind.  This  follows  a  rapid  change  in  the 
direction  and  magnitude  of  ion  convection,  which  is  not  immediately  reflected  in  the  neutral  flow 
because  the  neutral  atmosphere  responds  slowly  to  changes  in  the  ion  flow. 

This  result  shows  that  the  effect  of  neutral  winds  on  the  local  Joule  heating  rate  is 
nonnegligible  and  can  at  times  dominate.  During  most  of  the  time  illustrated,  the  calculated  Joule 
heating,  on  the  average,  is  2  to  4  times  smaller  when  the  neutral  wind  is  included  than  when  it  is 
set  to  zero. 

The  neutral  winds  have  seasonal  variations  (Roble  1983;  Wickwar  et  al.,  1986)  that  have  not 
been  considered  in  this  study.  These  variations  must  also  be  reflected  in  the  seasonal  variation  of 
Joule  heating. 
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SUMMARY  OF  OBSERVATIONS 

In  this  paper  we  have  focused  on  the  three  physical  parameters  -  the  electric  field,  the 

Pedersen  conductivity,  and  the  neutral  wind  -  that  determine  the  Joule  heating  rate.  Summarized 

below  are  some  significant  findings. 

•  The  data  presented  here  and  in  de  la  Beaujardtere  et  al.  (1988)  indicate  that  the  intensity  of 
the  ionospheric  electric  field,  as  well  as  the  average  convection  pattern,  are  different  for  each 
season.  The  magnitude  of  the  plasma  convection  is  smallest  in  summer  and  largest  in  fall. 

•  During  quiet  geomagnetic  conditions,  the  Joule  heating  rate  is  maximum  in  the  fall. 

•  The  F-region  contributes  little  to  nighttime  integrated  Pedersen  conductivities.  During  the 
daytime,  under  solar  maximum  conditions,  the  F-region  conductivity  can  be  larger  than  the  E- 
region  conductivity. 

•  The  often-neglected  neutral  wind  can  be  important.  In  our  example,  the  Joule  heating, 
calculated  by  incorporating  the  actual  neutral  wind  values,  was  generally  2  to  4  times  smaller 
than  that  calculated  by  assuming  U  =  0.  The  reverse  can  also  be  true,  as  was  illustrated 
during  a  period  when  the  wind  was  larger  than  the  E  X  B  drift. 


DISCUSSION 

An  inherent  limitation  in  averaging  all  data  by  season  as  we  have  done  here  is  that  we  smooth 
out  all  levels  of  magnetic  activity  and  of  IMF  Bz  and  By,  all  of  which  profoundly  affect  ionospheric 
plasma  convection.  However,  the  present  study  gave  enough  points  in  each  bin  so  that  the  results 
are  statistically  significant.  With  time  we  will  accumulate  enough  measurements  to  calculate 
separate  empirical  models  as  a  function  of  season,  magnetic  activity  and  IMF. 

As  shown  by  Equation  [2]  and  Figure  7b,  the  F  region  Pedersen  conductivity  depends  on  the 
collision  frequency  between  neutrals  and  0+  ions.  However,  it  is  becoming  increasingly  clear  that 
the  values  of  0+/0  collisional  cross  section  that  are  generally  adopted  may  have  to  be  altered.  As 
is  pointed  out  by  Burnside  et  al.  (1987)  the  value  commonly  adopted  may  have  to  be  multiplied  by 
as  much  as  1.7.  If  this  were  the  case,  the  F-region  Pedersen  conductivity  would  also  have  to  be 
increased,  since  0+  ions  dominate  at  this  altitude  (see  Figure  7a).  This  will  increase  even  more 
the  relative  contribution  from  the  F  region  to  the  integrated  Pedersen  conductivity. 

The  large  F  region  conductivity  has  important  consequences  in  terms  of  high  latitude 
electrodynamics:  (1)  Over  a  significant  fraction  of  the  F  layer,  the  ions  do  not  move  in  the  ExB 
direction,  but  rather,  collisions  have  the  effect  of  deflecting  them  in  the  direction  of  E.  (2)  The 
ionospheric  current  vector,  Jj_,  rotates  by  90°  between  the  E  and  F  layers.  Jx  is  perpendicular  to 
the  electric  field  in  the  E  region,  and  parallel  in  the  F  region.  (3)  A  large  fraction  of  the  field- 
aligned  current  closure  by  Pedersen  current  occurs  in  the  F  region. 

Our  last  point  concerns  the  role  of  the  neutral  wind  in  Joule  dissipation.  The  energy  dissipated 
to  the  neutrals  is  governed  by  Equation  [1],  and,  as  we  have  shown,  the  neutral  wind  term  should 
not  be  neglected  in  the  calculation.  When  dealing  with  the  magnetosphere,  however,  the  actual 
energy  dissipated  by  Joule  heating  is  evaluated  by  considering  only  the  electric  field  and  by  not 
including  the  neutral  wind  term  (Sugiura,  1984;  Vasyliunas,  1970). 
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Fig.  1.  Average  ion  convection  velocities  for  (a)  spring;  (b)  summer;  (c)  fall;  (d)  winter.  The 
Sondrestrom  data  from  March  1983  to  July  1988  were  averaged  into  1°  by  lhr  bins.  Local  time  is 
indicated  on  the  outside  and  UT  on  the  inside  of  the  clock  dial  circle. 


Fig.  2.  Components  of  the  average  ion  velocities  at  75°  invariant  latitude  with  (a)  eastward  and  (b) 
northward  velocities.  Note  that  the  velocity  scale  for  (a)  is  half  that  for  (b).  The  error  bar  indicates 
the  mean  standard  deviation.  Fall  and  equinox  have  been  averaged  together. 

Fig.  3.  Latitude  of  the  ion  convection  reversal  for  three  seasons.  This  reversal  is  that  of  the  2°X  2hr 
averaged  data.  Reversals  in  the  morning  and  afternoon  convection  cell  are  indicated  by  separate 
curves,  and  the  discontinuities  in  the  noon  and  midnight  local  times  show  where  the  transition 
takes  place.  Arrows  indicate  where  the  convection  reversal  is  located  beyond  the  radar  field  of 
view. 


Fig.  4.  Average  of  squared  ion  velocity  vectors  as  a  function  of  local  time  for  (a)  fall  and  summer 
(solid  and  dotted  lines,  respectively);  and  (b)  winter  and  spring  (thick  and  thin  lines,  respectively). 


Fig.  5.  Average  patterns  of  the  squared  electric  field  magnitude  measured  by  binning  AE  C  data 
(from  Foster  et  al.,  1983). 


Fig.  6.  Pedersen  height-integrated  conductivity  versus  time  measured  during  (a)  solar  minimum 
at  Sondrestrom  and  (b)  and  (c)  solar  maximum  at  Chatanika.  The  solid  curve  shows  the  Pedersen 
conductivity  integrated  up  to  a  500-km  altitude,  and  the  dashed  curve  is  for  a  180-km  altitude 
limit.  The  solar  zenith  angle  is  indicated  at  the  top  of  the  plots. 


Fig.  7.  Height  profiles  of  the  four  components  of  (a)  Hall  conductivity  and  (b)  Pedersen 
conductivity  (see  Equations  [2]  and  [3].)  The  data  correspond  to  14  LT  (00  UT)  on  February  4, 1981 , 
illustrated  in  Figure  6b. 


Fig.  8.  Electron  density  profile  corresponding  to  the  Figure  7  conductivities. 

Fig.  9.  Joule  heating  rates  at  a  115-km  altitude  calculated  neglecting  (solid  curve)  and  including 
(dashed  curve)  the  effect  of  the  neutral  wind.  The  data  correspond  to  June  11, 1980,  illustrated  in 
Figure  6c. 
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FIGURE  9.  JOULE  HEATING  RATES  AT  115  KM  CALCULATED  FOR  CHATANIKA  DATA  ON  800611  NEGLECTING 
(soM  curve)  AND  INCLUDING  (dashed  curve)  THE  EFFECT  OF  THE  NEUTRAL  WIND. 


